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We study a three-state {leftist, rightist, centrist) model that couples the dynamics of social balance 
with an external deradicalizing field. The mean-field analysis shows that there exists a critical value 
of the external field Pc such that for a weak external field {p<pc), the system exhibits a metastable 
fixed point and a saddle point in addition to a stable fixed point. However, if the strength of the 
external field is sufficiently large {p>pc), there is only one (stable) fixed point which corresponds 
to an all-centrist consensus state (absorbing state). In the weak-field regime, the convergence time 
to the absorbing state is evaluated using the quasi-stationary distribution and is found to be in 
agreement with the results obtained by numerical simulations. 
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I. INTRODUCTION 

Structural balance is considered to be one of the key 
driving mechanisms of social dynamics and since 

its introduction by Heider [I|, it has been studied ex¬ 
tensively in the context of social networks In a 

socially interacting population, relationships among in¬ 
dividuals {links in the underlying social network) can be 
classified as friendly (-1-) or unfriendly (-). Evolution of 
these links is governed by the theory of structural bal¬ 
ance, also referred to as social balance. The underlying 
axioms behind this theory are: (i) a friend of my friend 
or an enemy of my enemy is my friend, and (ii) a friend 
of my enemy or an enemy of my friend is my enemy. In 
the context of social networks, a triangle is said to be un¬ 
balanced if it contains an odd number of unfriendly links 
According to the theory of social balance, these un¬ 
balanced trian gle s have a tendency to evolve to balanced 
configurations [l|, 1^ . This might happen by transition¬ 
ing nodes and/or interpersonal links in such a way that 
the conditions of social balance are satisfied. For exam¬ 
ple, a triad in which two mutually antagonistic individ¬ 
uals have a common friend, is by definition unbalanced 
but can become balanced by requiring either the com¬ 
mon friend to choose a side or the mutually antagonistic 
two nodes to reconcile their conflict and become friends. 
A structurally balanced network contains no unbalanced 
triangles. Here, we construct an individual-based model 
where the dynamics (in part) is driven by structural bal¬ 
ance, but (unlike in previous works [5l4l0l|) a change in 
the state of an edge is the direct consequence of a change 
in the state (opinion) of one of the nodes the edge con¬ 
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nects. Another key feature of the model studied here is 
that triadic (three-body) interactions among nodes have 
been considered as opposed to dyadic (pair-wise) inter¬ 
actions in three-state models |11 | - Ih^ ). 


In this paper, we consider a population where each 
individual is in one of the three possible opinion states 
{leftist, rightist, or centrist) |lll-ll3l|. A link that connects 
two extremists of opposite type (i. e. the link between a 
leftist and a rightist) is considered to be unfriendly while 
all other links are friendly. Thus, a triangle containing 
one node of each type (leftist, rightist, and centrist) is 
unbalanced. An unbalanced triangle can be balanced in 
a number of ways with each minimal change solution re¬ 
quiring one node in the triangle updating its opinion. In 
a model with extremist and moderate opinion states of 
individuals. Marvel et al. [3 showed that moderation 
by external stimulus is a way to have a society adopt a 
moderate viewpoint (non-social deradicalization). Here, 
in our model, we consider a similar external influence field 
(e.g., campaigns, advertisements) which converts extrem¬ 
ists into centrists. Thus, the system is governed by the 
competing effects of social balance and influence. 


At each time step either with probability p (i) random 
node is selected and if it is an extremist, it is converted to 
a centrist, or with the complementary probability (1 — p) 
(ii) a random triangle is selected and if unbalanced, it 
is balanced by either converting (with a probability a) a 
centrist into an extremist or [with a probability (1—a)] an 
extremist into a centrist. Furthermore, since an extremist 
can either be a leftist or a rightist, a choice of converted 
extremist flavor is made with equal probability (i) as 
shown in Fig. [T] 
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FIG. 1. The balance step is taken with probability (1 — p). 
However, balance can be achieved in two ways: (a) with prob¬ 
ability a a centrist is converted into an extremist and since 
it can either be a leftist or a rightist, its flavor is chosen ran¬ 
domly with equal probability 1/2. (b) Alternatively, with 

probability (1 — tr) an extremist (either leftist or rightist with 
equal probability 1/2) is converted to a centrist. 


II. FULLY-CONNECTED NETWORKS 
(MEAN-FIELD ANALYSIS) 

A. Fixed points of the system 

For a fully-connected network, at any given time the 
state of a system of size N can be described by two num¬ 
bers - the density (fraction) of leftists (x) and the density 
of rightists (y) - as we can eliminate the density of cen¬ 
trists {z) since x + y + z = 1. Thus, the evolution can 
be mapped onto the xy plane. A finite system will al¬ 
ways have only one absorbing fixed point (for p > 0) 
that is a consensus state where every node has adopted 
the centrist opinion. First, for simplicity, we consider 
this dynamics on an infinite complete graph where every 
node is connected to every other node (i.e., in the mean- 
field limit). Starting from an arbitrary state (x > 0, 
p > 0, z > 0), in the absence of an external influenc¬ 
ing field (p = 0), the final state of the system is either 
polarized {z = 0) or is a coalition where mixed popula¬ 
tion of centrists and extremists of one kind (either leftist 
or rightist) coexists. For p = 0, the whole triangular 
boundary of the phase space in the xy plane becomes 
absorbing, thus, a pure consensus state in this case can 
not be reached through transitions from a different initial 
state (see Appendix A). For p > 0, the evolution of x,y 
densities is governed by the following rate equations: 

dx 

— =-px + 5{2a-l){l-p)xy{l-x-y) (1) 

^ =-ra + 3(2a - 1)(1-p)xp(l -X - y) . (2) 

A trivial solution of these equations is the all-centrist 
consensus state, i.e., (x,y) = (0,0) (or equivalently 2 = 
1). However, the steady state solution (see Appendix B) 
of these equations with a > ^ shows the existence of a 
critical point. 


3(2a- 1) 
8-b3(2a- 1) ’ 


(3) 


such that for p < pc the system exhibits two non-trivial 
fixed points as well: 

which is a metastable fixed point and 

which is a saddle point (unstable fixed point). In the 
other scenario, when a < |, the system already has the 
tendency to move towards an all-centrist consensus state, 
hence only the trivial fixed point (x, y) = (0, 0) exists. 

In this paper we focus on the case of a > ^ in which 
the system has the tendency to become polarized and the 
external influence field is required to prevent this polar¬ 
ization. In this regime, due to the competition between 
balancing and influencing forces, the densities fluctuate 
around the metastable point and the system is trapped 
for exponentially long times before unlikely large fluc¬ 
tuation moves it to the absorbing state. Here all fixed 
points lie on the line y = x as shown in Fig. [2] and any 
asymmetry in x and y decays exponentially fast [l5|, [l6| 
(see Appendix B). The trajectories shown in Fig. are 
exact only in the thermodynamic limit (fully-connected 
network in the limit of Y —>■ oo). It can be seen explicitly 
that a finite network {N = 100) in the weak-held limit 
(p < Pc) gets stuck in the metastable state and never 
crosses the saddle point whereas a fast centrist consen¬ 
sus is reached when p > Pc (Fig. [3]). It is also clear from 
the above that the locations of the fixed points (roots) 
depend on the choice of a (> ^) and p. For a particular 
choice of a, the two fixed points (metastable and saddle) 
move closer to each other as p is increased from 0 till 
they meet and annihilate each other at p = pc (as shown 
in Fig.|31). Beyond Pc, these additional fixed points cease 
to exist and the only fixed point is the consensus state 
[x = y = 0). 

B. Consensus time for finite-size networks 

An all centrist consensus state is always reached for 
a finite network. Time to reach this absorbing state 
(consensus time Tc) can be obtained by direct simula¬ 
tions. This approach works well for p > pc, however, for 
p < Pc (specially when p << pc and/or N >> I), Tc 
becomes so large that its estimation by simulation be¬ 
comes difficult if at all possible. We therefore use the 
quasi-stationary ( QS) approximation prescribed in da 
and also used in (iSl Il9j to estimate Tc in the region 
P <Pc- 

We start by introducing notation for the numbers of 
nodes with the given opinion, thus, X = xN, Y = yN, 
and Z = zN. Then we form the master equation that de¬ 
scribes the time evolution of probability Px,y (the prob- 
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FIG. 2. (Color online) Phase-space trajectories within mean- 
field approximation for Q!=0.75 for which pc~0.16, (a) for 
p=0.12<pc, and (b) for p —0.20>pc- 



X 



X 

FIG. 3. (Color online) Stochastic trajectories for a fully- 
connected network of size N — 100 and a = 0.75, (a) for 
p=0.12<pc, and (b) for p =0.20>pc. See also Snpplemental 
Material for animations of the above two scenarios (for the 
same parameters). 


ability that system has X leftists and Y rightists at time 
t). 



p 


FIG. 4. (Color online) The locations {x, y co-ordinates) of 
metastable (solid lines) and saddle (dashed lines) fixed points 
for a = 0.75 (red) and a = 1.0 (green) obtained from the 
solution of the rate equations, x = y for all fixed points. 
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Within the triangular region (bounded by 0 < X < 
{N — Y) and 0 < F < {N — X)), the transitions al¬ 
lowed from a state (X, Y) are to states {X ± 1, Y) or 
{X, Y ±1) with the constraint that the system stays 
within the bounded region. The positive and negative 
terms on the right side of the master equation contribute 
to the net flow of probability into and out of the state 
{X,Y), respectively. A factor of ^ on the left-hand side 
appears because a microscopic step of transition from ini¬ 
tial state to the final state takes place in a time interval 
1/A^. ^ N(N density of unbalanced trian¬ 
gles at any given time for a fully-connected network. 

The QS distribution of occupation probabilities is 
given by Px,y = Px,Y{t)/Ps{t) where Ps{t) is the sur¬ 
vival probability. Under the QS hypothesis, the survival 
probability decays exponentially, governed by 

= -Psmo , ( 7 ) 

where Qq = _p[Ui,o + Po,i] measures the flow of probabil¬ 
ity into the absorbing state (0,0). The underlying idea 
of QS hypothesis is that the occupation probability dis¬ 
tribution conditioned on survival (over all {X,Y) except 
the absorbing state) is stationary. Therefore, 


convergence was obtained in 40000 iterations. As ex¬ 
pected from the mean-field analysis (for a = 0.75, the 
metastable fixed point (a:, y) = (0.38, 0.38)) the distribu¬ 
tion peaks around (X, Y) = (38, 38). 



100 

X 


FIG. 5. (Color online) The QS distribution with a=0.5 for 
A=100, a=0.75, p=0.12. For these parameters, pc~0.16. 


dPx,Y 5 dPs{t) 


( 8 ) 


We plug in Px,y in terms of Px,y into the master 
equation to obtain the QS distribution, 


Qx,y 

Wxx — Qo 


(9) 


where Wx,y = 

3a(l -p){X - 1)F(X -X + l-Y) 


Qx,Y = Px-l.Y 
+ Px,Y-l 
+ Px+l.Y 


{N -1){N -2) 

3a(l - p)X{Y -1){N - X - Y + 1) 

{N -1){N -2) 

3(1 - a)(l - p){X + 1)F(X -X -l-Y) 


(N -1){N -2) 

+ Px+i.y P {X + 1) 

3(1 - a)(l - p)X{Y +1){N-X-Y-1) 


+ Px, 


Y+l- 


{N-l){N-2) 


+ Px,Y + l P {Y + 1) . 


Starting from an arbitrary distribution P^ y i nn 
asymptotic QS distribution Px,y can be obtained by the 
iteration; P^y = a Pxy + (1 ~ n) i where 

’ ’ V Qq 

0 < a < 1 is an arbitrary parameter [2. The QS dis¬ 
tribution for a particular system size N = 100 (fully- 
connected), p = 0.12, and a = 0.75 with parameter 
a = 0.5 is shown in Fig. [5] In this case a satisfactory 


Once the desired QS distribution is obtained, the mean 
consensus time Tc is computed from the decay rate of the 
survival probability. 




P [Pi,o + Po,i] 


( 10 ) 



FIG. 6. (Golor online) Gonsensus time Tc obtained from the 
QS distribution and by direct MC simulations for q=0.75 

(pc« 0 . 16 ). 

We compare Tc obtained from the QS approximation 
to that obtained by direct Monte Carlo (MC) simulations 
in the region of p < Pc where Tc could be easily obtained 
by both methods for the respective system sizes [Fig.|6]. 
One can see that there is a good agreement between the 
two methods across many system sizes and the agreement 
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is getting better as p is decreased below pc as shown in 
Fig. [HI The change in the relative error with respect to 
p and N, e = I ^ can be seen in the inset. 
Figure [7] shows Tc obtained by the QS approximation as 
a function of p for the entire range of p considered for all 
system sizes, where MC simulations become prohibitive 
to estimate T^. The consensus time Tc shows an expo¬ 
nential scaling with N for p < pc as shown in Fig. |S1 



FIG. 7. (Color online) Consensus time Tc computed by the 
QS approximation as a function of p in the weak-field regime 
for q=0.75 (pc«0.16). 



FIG. 8. (Color online) Consensus time Tc computed by the 
QS approximation as a function of system size N in the weak- 
field regime for a=0.75 (pc«0.16). 

In order to obtain the dependence of Tc on p, we as¬ 
sume a relation common in systems with tipping points 
and barrier crossing BUI, 

Tc = f{N)exp\f3{p)N], (II) 

where f{N) is increasing slower than exponential with N 
and 

/5(p) ~ be-^r • (12) 

With Tc{pc) = /(IV), we have 

In(Tc) - ln(Tc(pc)) = fi{p)N . (13) 


As can be seen from Fig. |9l the growth rate approxi¬ 
mately follows the scaling behavior /3 ~ \Pc — pV with 
the measured exponent v « 1.57. 



FIG. 9. (Color online) In(rc) — ln(rc(Pc)) as a function of 
(pc—p) on a log-log scale. The exponent v is given by the slope 
of the fitted line {N = 300, a = 0.75, pcSsO.16, v « 1.57). 


III. LOW-DIMENSIONAL NETWORKS 

For sparse low-dimensional networks, the mean-field 
analysis does not hold and QS approach is difficult to for¬ 
mulate. Hence we rely solely on MC simulations. Specif¬ 
ically, we look at the survival probability Pg (for a hxed 
cutoff time t = 5000) for a 2D random geometric graph 
(RGG) [m with (fc) = 10 and a ID regular lattice with 
each node having a degree k = 10. We start with a 
polarized initial state (x = 0.5,?/ = 0.5). The simula¬ 
tion results indicate the existence of a critical point pc at 
which the survival probability undergoes an abrupt tran¬ 
sition [ 2 ^ as shown in Fig.[T0l We choose these particular 
spatial embeddings because presence of local clustering 
ensures a significant number of triangles in the network 
and the model requires the presence of triangles for the 
balance dynamics to take place. For networks with rela¬ 
tively low clustering coefficient (e.g. ER, BA networks), 
the dynamics would be heavily dominated by external 
influence. To examine the dynamics on a real-world net¬ 
work, we also simulated the dynamics, starting from the 
same initial conditions (randomly assigning the opinions 
in the initial state of the system), on the giant component 
of a high-school friendship network from the Add Health 
data set . The critical point in pc is shown to exist in 
this network structure as well (Fig. ITO)) . 

The simulation results (Fig. fTO)) indicate that the criti¬ 
cal point for 2D RGG is significantly higher than that for 
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FIG. 10. (Color online) Survival probability at time t=5000, 
as a function of p for Q!=1.0 for a ID regular lattice (with 
degree fc=10), a 2D RGG (with (A:) = 10) with A'^=1000, and 
for a high-school friendship network with (fc)«6 and N = 921. 


ID lattice. For the same choice of a = 1, the critical point 
in the case of fully-connected network is pc « 0.27 (be¬ 
tween the values of the ID and 2D systems). Simulation 
results plotted in Figs. fTlT a) and (b) show that the decay 
of fraction of unbalanced triangles in the network (n„) is 
governed by the power-law in the case of 2D RGG but it 
decays exponentially (ignoring the transience) in ID reg¬ 
ular lattice as shown in Figs, fl^ a) and (b). In the case 
of 2D RGG, there are frustrated domains that are long- 
lived (shown in Fig. [T51) . whose existence can be causal 
or symptomatic to the slow convergence. Particular net¬ 
work structure and spatial correlations (among nodes) 
become significantly important in defining the dynamics 
of low dimensional systems and a quantitative analysis 
becomes mathematically challenging. 


IV. SUMMARY 

In summary, we presented a framework that models 
social imbalance arising from individual opinions in the 
simplest manner possible, and observe its counteracting 
effect on an externally influencing field. We found that 
there exists a critical value Pc, above which the only fixed 
point is the centrist consensus state, and below which a 
metastable fixed point of the system emerges. We demon¬ 
strated how the competition between balance and influ¬ 
ence can lead the system to metastability. Using a semi- 
analytical approach (QS approximation), we estimated 
the consensus times which show good agreement with 
simulation results. Additionally, employing simulations, 
we demonstrated that this critical behavior is also seen 
in sparse networks. 




FIG. 11. (Golor online) Average fraction of unbalanced tri¬ 
angles in the 2D RGG with N — 1000, (fc) = 10, and a = 1.0 
(a) on a log-log scale and (b) on a semi-log scale. The initial 
population densities are x = 0.5, y — 0.5. 
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Appendix A: The p = 0 case 

In this appendix we discuss the p = 0 case of the model 
(case with no external influence). When p = 0, the rate 
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FIG. 12. (Color online) Average fraction of unbalanced tri¬ 
angles in the ID regular lattice with N = 1000, k = 10, and 
a = 1.0 (a) on a log-log scale and (b) on a semi-log scale. The 
initial population densities are x = 0.5, y — 0.5. 


equations for the densities x and y can be written as 

(IT 

— = 3{2a - l)xy{l - x - y) 

^ = 3{2a - l)xy{l - X - y) . (Al) 

at 

In this case, the entire boundary of the triangular phase 
space becomes absorbing because structural balance is 
achieved as soon as either of the x, y, or z variables 
becomes zero and from that point the system does not 
evolve. The type of steady-state (in other other words, 
which of the three boundaries is hit by the system) of 
this system depends on the choice of a. For a < ^, 
the system moves towards an all-centrist consensus till it 
eventually hits either a; = 0 (centrist-leftist coalition) or 
y = 0 (centrist-rightist coalition) boundary. However, for 
a > 5 (in the absence of external field, p = 0) the system 
tends to more radical configuration and stops evolving 
when x + y = 1 (or reaches the long side of the triangular 
phase space). After gaining some insight in this case, 
we performed stochastic simulations. The probability for 
the system to end up in a polarized state Plr is shown 
in Fig. [T3] (with the complementary probability Pce = 
1 — Plr, the system settles in a coalition state). Starting 



FIG. 13. (Color online) Time of evolution (single-run) of the 
system. Green, red, and blue nodes correspond to leftists, 
rightists, and centrists respectively. The network structure is 
2D RGG with (fc) = 10, = 1000, a = 1.0, and p = 0.3. 
System is initialized with x — 0.5, y = 0.5. (a) t = 0, (b) 
t = 10, (c) t = 100, and (d) t = 500, where t is the number of 
time steps. 


from an equal density initial state (xq ~ yo ~ ~ 0.33), 

we also look at the composition of the final state when 
the system reaches a polarized [a = 0.8) or a coalition 
state (a = 0.3). We only show the distribution of x 
given that the system reaches a final state on the y = 













Appendix B: Steady-state solution of the rate 
equations 

The rate equations for the leftist and rightist densities 
under the mean-field assumption are 

dx 

— = -px + 3(2q; - 1)(1 - p)xy{l - x - y) 

^ = -py + 3(2a - 1)(1 - p)xy{l - x - y) . (Bl) 
at 

By adding and subtracting the above equations and intro¬ 
ducing a new set of variables u = {x + y) and v = (x — y), 
one can immediately see that 


FIG. 14. (Color online) The probability Plr as a function of 
a for N = 100 (fully-connected) and different starting points 

(p = 0). 


dv 

dt 


= -pv , 


(B2) 


yielding v ^ exp{—pt), which means that ?; —>■ 0 exponen¬ 
tially fast. Therefore, we can assume that x ~ y which 
allows us to analyze the system in terms of single-variable 
equation for the “slow” mode u (and x = y = m/2), 


0 or 2 / = 1 — X boundary (Fig. [15)). y has the same 
distribution along x = 0 and y = 1 — x boundaries due to 
symmetry. Starting from an arbitrary state (xq > 0, j/o > 
0, zq > 0), the system never reaches a pure consensus 
because a structurally balanced configuration is always 
reached before reaching a pure consensus state and the 
system freezes in that state. The consensus in this case 
is observed only if the initial state itself is a consensus 
state and the system remains in that state. 


du 

dt 


—pu + 6(2a 


1)(1 - P ) 




(B3) 


which can be solved for stead-state ^ = 0. A trivial 

at 

solution of this equation is m = 0, which is the absorbing 
state (x = 0, y = 0). Additional roots are the solutions 
of the quadratic equation 


M^ — M -I- 


2p 

3(l-p)(2a 



(B4) 


and are given by 



X (leftist density) 


= i + i /. 

^ 2 2\j 3(1-p)(2a-1) ■ 


(B5) 


These solutions make sense only when a > 1/2, otherwise 
the solution will lie outside the feasible domain [(x-|-y) < 
1]. For a > 1/2, we obtain a critical point, 


3(2a- 1) 
8-b3(2a- 1) ’ 


(B6) 


such that the roots are real and positive for p < pc- Thus, 
in terms of x and y the two roots (other than the absorb¬ 
ing state) are 


FIG. 15. (Color online) The histogram of leftist density x 
when the system ends up in a leftist-rightist polarization or 
centrist-leftist coalition state for a = 0.8 and a = 0.3. The 
network is fully-connected {N = 100 and p = 0). 


^ ^ 3(l-p)(2a-l) 

1 1 

^ ^ 4 4y 3(l-p)(2a-l) 


(metastable) 


(saddle) .(B7) 
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